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Abstract 

Computer algebra programs are presented for application in general relativity, in 
electrodynamics, and in gauge theories of gravity. The mathematical formalism used 
is the calculus of exterior differential forms, the computer algebra system applied 
Hearn's Reduce with Schriifer's exterior form package Excalc. As a non-trivial ex- 
ample we discuss a metric of Plebahski & Demiahski (of Petrov type D) together 
with an electromagnetic potential and a triplet of post-Riemannian one- forms. This 
whole geometrical construct represents an exact solution of a metric-affine gauge 
theory of gravity. We describe a sample session and verify by computer that this 
exact solution fulfills the appropriate field equations.- Computer programs are de- 
scribed for the irreducible decomposition of (non-Riemannian) curvature, torsion, 
and nonmetricity. file cpc9.tex, 1998-03-18 



1 Introduction 



General relativity (GR), that is, Einstein's gravitational theory, is one of 
those scientific subjects in which computer algebra methods were used as 
soon as they became available. In fact, in some instances, general relativ- 
ity was taken as a guinea pig since this theory is notorious for its lengthy 
and messy calculations. Just remember that the Riemann curvature 2-form 
Ra/3 '■= 9/3-yRa 1 = —R/3a, here are the components of the metric, char- 
acterizes the local geometry of four- dimensional spacetime. It has, in these 
dimensions, 20 independent components. And Einstein's field equation has 
10 independent components. These numbers alone make it plausible that it 
is desirable to give the corresponding calculations to a computer, see Brans 
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[3]. Since they are analytic, numerics doesn't help, rather symbol manipulation 
programs or, as they are called today, computer algebra systems are necessary. 

This is even more important for generalizations of Einstein's theory, like the 
Einstein-Cartan theory, the Poincare gauge theory, or the metric-affine the- 
ory of gravity. In these gauge theories of gravitation, see Gronwald & Hehl 
[11] for a review, spacetime carries additional post-Riemannian structures, 
namely torsion T a (a vector valued 2-form) with 6x4 independent compo- 
nents, and/or nonmetricity Q a p = Qp a (a symmetric tensor- valued 1-form) 
with 4 x 10 components. Moreover, in general, the curvature loses its anti- 
symmetry. Thus, a metric-affine space, in four dimensions, is described by 
the 36 + 24 + 40 'field strengths' R a/ 3, T a , and Qap, respectively. For some 
physicist this is too much and considered to be a plethora of arbitrariness. 
This is like the argument against Einstein's theory, which was taken quite 
seriously in some quarters, that the 10 gravitational potentials of Einstein's 
theory are too much as compared to only 1 potential in Newton's gravitational 
theory. But consistency and beauty cannot be measured in numbers of com- 
ponents! In [16], e.g., one can find the arguments which make us believe that 
the microstructure of spacetime requires post-Riemannian degrees of freedom 
as represented by torsion and nonmetricity. 

Let us come back to general relativity. A readable survey of the computer 
algebra systems in use for general relativity has been given by Hartley [12]. 
In the same book [17], see also [14], there is a whole chapter on computer 
algebra methods as used in relativity. Suppose you would like to use such 
methods in general relativity or gauge theories of gravity. Considering the 
ample computer resources one has available today, it seems reasonable to 
us to turn to a general purpose system like Macsyma, Maple, Mathematica, 
and/or Reduce and, on top of one of these, to use specialized packages for 
general relativity/differential geometry, see [12]. For getting familiar with the 
corresponding general relativity packages, we can recommend the respective 
articles or books of McLenaghan [23] for Maple, of Parker & Christensen 
[26] and Soleng [34] for Mathematica, and of McCrea [22] for Reduce. Some 
specific relativistic applications have been discussed, e.g., in Maple, see [28], 
in Mathematica, see [19,39], and in Reduce, see [43]. 

The general relativity computer algebra packages can, as a rule, also be used 
for non-Riemannian spacetimes and the corresponding gravitational gauge the- 
ories. Usually the inclusion of torsion and nonmetricity is not difficult to pro- 
gram. However, there is not too much literature available on that subject. For 
torsion theories there exist, amongst others, the papers of Aman et al. [1,2], 
McCrea [22] , Schriifer et al. [33] , Tertychniy and Obukhova [36] , and Zhytnikov 
[44]. The packages GRG [44] and GRG EC [36] (EC stands for Einstein-Cartan 
theory) are Reduce based systems. 
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There is even less literature available for metric- a ffine spaces and the cor- 
responding gauge theories. Such a spacetime carries, besides the coframe (of 
1-forms) d a = ei a dx l , the metric g = g a p q d OL ®'d 13 = gij dx l ®dx : ' and the linear 
connection 1-form Yj 3 = T^a 13 -# 7 = Y io P dx % as totally independent gravita- 
tional 'potentials': (g a p, $ a , Fa 13 ). For the metric-affine gauge theory of gravity 
(MAG), which encompasses the Poincare gauge theory, the Einstein-Cartan 
theory, and general relativity as specific subcases, we developed a couple of 
computer algebra programs which were and are used in our respective groups 
for the search of exact solutions, e.g. We would like to explain and demonstrate 
them. They were developed over the last eight years or so and different people 
were involved: Dermott McCrea^, Werner Esser, Frank Gronwald, Ralf Hecht, 
Yuri Obukhov, Roland Puntigam, Sergei Tertychniy, Romualdo Tresguerres, 
Eugen Vlachynsky, and, more recently, Jose Socorro. 

We will proceed as follows: The general metric-affine framework, including all 
the conventions, will be taken from [16], see also [11]. Our computer algebra 
programs use Cartan's calculus of exterior differential forms as implemented 
on the computer by means of Hearn's computer algebra system Reduce [13] 
and Schrufer's corresponding Excalc package [32]. We will always have in mind 
the special case of general relativity, i.e., we can specialize to general relativ- 
ity if desirable. This is important if one wants to extract the pos^-Riemannian 
structures from the computed expressions, i.e., the deviations from the corre- 
sponding Riemannian spacetime. We are trying to program our MAG proce- 
dures as some kind of master programs which can be immediately and effi- 
ciently specialized to Riemann-Cartan or Riemann spacetimes. 

In this sense, we turn, in Sec. 2, to the curvature 2-form of a metric-affine 
space with its 16 x 16 independent components. We will irreducibly decom- 
pose it under the Lorentz group and discover thereby typical post-Riemannian 
pieces, besides the pieces which emerge already in a Riemannian space. This 
automatized irreducible decomposition of the curvature may also be useful for 
people who are not concerned with gravitational gauge theories but just with 
non- Riemannian structures in differential geometry, i.e., we believe that Sec. 2 
should be of more general interest than the rest of the paper. 

In Sec. 3, as an example for testing our methods, we turn to an exact (elec- 
trovac) solution of the Einstein-Maxwell equations of Petrov type D as formu- 
lated by Plebanski and Demianski (P&D) [27], see also Garcia Diaz [9] and 
Debever, Kamran, and McLenaghan [4,5]. After a short review of the Einstein- 
Maxwell equations in exterior calculus, we develop the computer code for the 
orthonormal coframe of the P&D solution and compute its curvature. 

Subsequently, we will determine the electromagnetic energy-momentum distri- 
bution of the P&D solution. For this purpose we bring in our Maxwell program 
code which is suitable for electromagnetic fields on top of an arbitrary differ- 
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entiable manifold, whatever metric and whatever connection, see [30,29]. This 
so-called metric-free formulation of Maxwell's equations (the metric enters 
only via the Hodge star into the constitutive assumption) is extremely concise 
and well adapted to executing such computations within geometrical theories 
of gravitation. 

We combine the decomposition method of Sec. 2 with this Maxwell program. 
We can now easily verify the Einstein-Maxwell equations and find, in addition, 
a very compact form for Weyl's conformal curvature 2-form. 

In Sec. 4, we turn our attention again to non-Riemannian geometries. We de- 
scribe programs for the irreducible decomposition of torsion and nonmetricity. 

In Sec. 5, the P&D metric of GR is extended to an exact solution in the frame- 
work of metric- affine gravity. For that purpose a non-Riemannian covector 
triplet has to be introduced, see [37,38,41,25,42,6,29,24,40,15]. Using our de- 
composition programs as described in Sees. 2 and 4, we are able to test this 
whole geometro-physical construct provided we also use our Maxwell package. 
It should be pointed out that this exact solution of metric-afline gravity was 
only found by means of an extensive use of computer algebra. In Sec. 6 we will 
conclude with a discussion. 



2 Curvature in non-Riemannian and Riemannian spaces, its irre- 
ducible decomposition by machine 

Before we turn to the programs, let us report on the irreducible decomposition 
of the curvature 2-form 



see [16] for a discussion and the relevant literature. Besides the linear con- 
nection r^, we have a metric g available with signature ( — h ++). We may 
lower the second index on the curvature two-form. We consider the irreducible 
decomposition of R a $ in 4 dimensions under the Lorentz group. The first step 
is to separate it into its antisymmetric and symmetric parts, 



R a/ 3 = W aP + Z a(3 ~ curv2(-a, -b) := w2(-a, -b) + z2(-a, -b) , (2) 



W a p := R [aPi ~ w2(-a, -b) := (curv2(-a, -b) - curv2(-b, -a))/2 (3) 




with 
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and 



Z a p := R[ a/3 ] ~ z2(-a, -b) := (curv2(-a, -b) + curv2(-b, -a))/2 . (4) 



We displayed at the same time the computer names (identifiers) of the quanti- 
ties emerging and the formulae in Reduce- Excalc syntax. A tutorial on the use 
of Excalc is contained in [35], e.g. We recall that lower indices are denoted by 
a minus sign. For convenience, in Excalc, we use Latin indices for the frames 
instead of Greek ones as in the formulae. The rank p of a form (if p > 0) is 
usually indicated by a corresponding Arabic number in the computer name 
of the quantity, that is, the curvature is a 2-form, thus we call it curv2. The 
Einstein form is a 3-form or einstein3(a), etc. 



2.1 Rotational curvature 



The antisymmetric piece W a p, the 'rotational curvature', can be decomposed, 
similarly as in a Riemann-Cartan spacetime, into 6 pieces, 

w aP = ww aP + {2) W a p + ^w a(3 + ^w a(3 + ^W a p + ^W a p 

~ weyl + paircom + pscalar + ricsymf + ricanti + scalar . (5) 
Here weyl corresponds to (%„^ etc. 

We can 'contract' the 2-form W a p in different ways by means of the frame e a ~ 
e(— a), the coframe ti a ~ o(a), and the Hodge star * ~ #. The corresponding 
truncated quantities read: 



W a := ep\W ap ~ wonel(a) 

W := e a \W a ~ wzero 

X a := *(W fhx A tip) ~ xonel(a) 

X := e a \X a ~ xzero 

They can further be decomposed into 



= e(-b)_|w2(a,b), (6) 

= e(— a)_ |wonel(a) , (7) 

= #(w2(b,a) Ao(-b)), (8) 

= e(— a)_ |xonel(a) . (9) 



V a := X a -±0 a AX- l -e a \ A Xp) , 
® a ■= W a - ^W$ a - X -e a \ A Wp) , 



(10) 
(11) 



or 



psil(-a) := xonel (-a)-(l/4) *o(-a) "xzero 
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-(l/2)*e(-a)_| (o(b)"xonel(-b))$ 
phil(-a) := wonel(-a)-(l/4)*wzero*o(-a) 

-(l/2)*e(-a)_| (o(b)"wonel(-b))$ 

and then 'blown up' to the 2-forms 



paircom2(-a, -b) *(# [q A tffl) , (12) 

pscalar2(-a, -b) ~ —^*(X A$ a A&p) , (13) 
ricsymf 2(-a, -b) % A <S> p] , (14) 

ricanti2(-a, -b) ~ ~-& [a A e p] \ (f> A W 7 ) , (15) 
scalar2(-a, -b) ~ — -Wti a A dp , (16) 

weyl2(-a, -b) ~ W afi - Y}^W afi • (17) 

1=2 

To be concrete, let us display (12), e.g., as a Reduce-Excalc command: 

paircom2(-a, -b) : =-(1/2) *#(o (-a) "psil (-b) -o (-b) "psil (-a) ) ; 
The hash sign # represents the Hodge star operator. 

Analogously, we can also translate the remaining definitions of the irreducible 
components into Reduce-Excalc. We remember that the rank of a form in 
Reduce-Excalc has to be declared explicitly by means of the pf orm declaration. 
In this way we arrive at the decomposition file: 

% Irred. decomposition of the rotational curvature w2(-a,-b) * 
°/ file rotcurvdecomp . exi , 1998-03-15 

% input: coframe o(a), frame e(a), metric g(a,b), curv2(a,b) 

pform {wzero ,xzero}=0 , 

{wonel (a) ,xonel (a) , psil (a) ,phil (a) }=1 , 

{w2(a,b) ,paircom2(a,b) ,pscalar2(a,b) , ricsymf 2(a,b) , 

ricanti2(a,b) , scalar2(a,b) , weyl2(a,b)}=2$ 



w2(-a,-b) 
wonel (a) 
wzero 
xonel (a) 
xzero 
psil(-a) 



= (1/2) * (curv2 (-a, -b) -curv2 (-b , -a) ) $ 

= e(-b)_|w2(a,b)$ 

= e(-b)_ | wonel (b)$ 

= #(w2(b,a)~o(-b))$ 

= e(-a)_|xonel(a)$ 

= xonel(-a)-(l/4)*o(-a) "xzero 
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-(l/2)*e(-a)_| (o(b)~xonel(-b))$ 
phil(-a) := wonel(-a)-(l/4)*o(-a)*wzero 

-(l/2)*e(-a)_| (o(b)"wonel(-b))$ 



paircom2(-a, -b) 
pscalar2(-a, -b) 
ricsymf 2(-a, -b) 
ricanti2(-a, -b) 

scalar2(-a,-b) 
weyl2(-a,-b) 



:= -(l/2)*#(o(-a)~psil(-b)-o(-b)~psil(-a)) ; 
:= -(l/12)*#(xzero~o(-a)~o(-b)) ; 
: = - (1/2) * (o (-a) ~phil (-b) -o (-b) "phi 1 (-a) ) ; 
:= -(l/4)*(o(-a)~ (e(-b)_| (o(c) "wonel (-c) ) ) 

-o(-b)~(e(-a)_| (o(c)-wonel(-c)))) ; 
:= -(l/12)*wzero*o(-a)~o(-b) ; 
: = w2 (-a, -b) -paircom2 (-a, -b) -pscalar2 (-a, -b) - 
ricsymf 2 (-a, -b) -ricanti2(-a, -b)-scalar2(-a, -b) ; 



clear wzero,xzero,wonel(a) ,xonel(a) ,psil(a) ,phil(a)$ 
$end$ 



Before this program segment can be used, we have to specify a coframe o(a), 
the frame e(-a) dual to it, a metric g, and a connection gamma(-a,b), and 
to compute the curvature curv2(-a,b) beforehand. This will be done below. 

Already in a Riemann-Cartan space (with metric compatible connection), i.e., 
in the case of vanishing nonmetricity, there emerge these 6 pieces of the cur- 
vature tensor with their 10 + 9+1 + 9 + 6 + 1 independent components. In a 
Riemannian space, we have paircom = pscalar = ricanti = 0, or 

-Ra/3 7 5 = R-ySa/3 , R[a/3jS\ = , RiC[a/3] = , (18) 



with Ric a /3 := i? 7Q/ g 7 . Then only weyl, ricsymf (ricci symmetric and tracef ree), 
and scalar survive with their, as is well known from general relativity, 10 + 
9 + 1 independent components. In other words, the torsion T a , via the first 
Bianchi identity DT a = R p a A induces the non-vanishing of paircom, 
pscalar, and ricanti. Obviously, (18) represents the 16 algebraic identities 
of a Riemannian curvature tensor. 



2.2 Strain curvature = segmental curvature © shear curvature 



The symmetric part Z a/3 of the curvature two-form, the 'strain curvature', is 
more involved than the rotational curvature. First of all, we split it into a 
tracefree (or shear) and a trace (or dilation) part: 

Z a p = X*/3 + \ 9af3 Z , Z := ZJ . (19) 
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The dilation piece g a pZJ /A is Weyl's segmental curvature, which we recognize 
as one irreducible piece of the strain curvature. The rest, the shear curvature 
/& a p, will be 'contracted': 

K--=e"\K P , A:=i(^A^J, Y a := *(^A^). (20) 

An even finer decomposition is possible: 

2. := K ~ \e a \ A ?J , Y Q := Y a ^ e Q j (f A Y,) . (21) 

In terms of these newly introduced quantities, the irreducible decomposition 
of Z a p may be written as 

Z a p = {1) Z a(S + ^Z a/3 + ^Z aP + ^Z aP + ^Z a(S , (22) 
~ zcurvone + zcurvtwo + zcurvthree + dilcurv + zcurvf ive , 

where 



zcurvtwo2(— a, 


-b) 


~ -^Kat,)}, 


(23) 


zcurvthree2(— a, 


-b) 


~ i {2%A(e /3) JA)-( ?Q/3 A}, 


(24) 


dilcurv2(— a, 


-b) 


1 7 

~ ^ 9af3 Z , 


(25) 


zcurvf ive2(— a, 


-b) 




(26) 


zcurvone2(— a, 


-b) 


5 


(27) 



1=2 



The corresponding Reduce-Excalc code reads as follows: 

% Irred. decomposition of the strain curvature z2(-a,-b) * 
% file straincurvdecomp . exi , 1998-03-15 

% input: coframe o(a), frame e(a), metric g(a,b) , curv2(a,b) 

pf orm {ztracef 1 (a) ,yyl (a) ,xil (a) ,upsilonl (a) }=1 , 
{z2(a,b) , ztracef 2(a,b) ,delta2, 
zcurvone2(a,b) , zcurvtwo2(a,b) ,zcurvthree2(a,b) , 
dilcurv2(a,b) , zcurvf ive2 (a,b)}=2$ 
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z2(-a,-b) 
ztracef 2(-a,-b) 
ztracef l(-a) 
delta2 
yyl(-a) 
xil(-a) 
upsilonl (-a) 



: = ( 1/2) * (curv2 (-a , -b) +curv2 (-b , -a) ) $ 
:= z2(-a,-b)-(l/4)*g(-a,-b)*z2(-c,c)$ 
:= e(b)_|ztracef2(-a,-b)$ 
:= (l/2)*o(a) "ztracef l(-a)$ 
:= #(ztracef2(-a,-b)"o(b))$ 

= ztracef l(-a)-(l/2)*e(-a)_| (o(c) "ztracef 1 (-c) )$ 
:= yyl(-a)-(l/2)*e(-a)_| (o(c)"yyl(-c))$ 



zcurvtwo2(-a,-b) := - (1/4) *#(o (-a) "upsilonl (-b) 

+o (-b) "upsilonl (-a) ) ; 
zcurvthree2 (-a, -b) : = (1/6) * (2* (o (-a) ~ (e (-b) _ I delta2) 

+o (-b) ~ (e (-a) _ | delta2) ) -2*g (-a, -b) *delta2) ; 
= (l/4)*g(-a,-b)*z2(-c,c); 
= (l/4)*(o(-a)"xil(-b)+o(-b)"xil(-a)) ; 
= z2(-a,-b)-zcurvtwo2(-a,-b) 



dilcurv2(-a,-b) 
zcurvf ive2(-a, -b) 
zcurvone2 (-a , -b) 



-zcurvthree2 (-a, -b) -dilcurv2 (-a, -b) -zcurvf ive2 (-a, -b) ; 



clear ztracef 1(a) , yyl(a) ,xil(a) , upsilonl (a) , 
ztracef 2(a,b) ,delta2$ 

$end$ 



3 Intermezzo: Einstein-Maxwell equations, Plebanski & Demianski 
metric 



Let us now turn to general relativity (GR) in order to have a decent electrovac 
solution of GR as a starting point for our further considerations, i.e. an ex- 
act solution of the Einstein equation with a distribution of electromagnetic 
energy-momentum as source. The Einstein-Maxwell equations with cosmolog- 
ical constant A read: 

Max 

G a + X Va = 2 T a , dH = , dF = 0, H = *F. (28) 

The Einstein 3-form G a ~ einstein3(— a) and the electromagnetic energy- 
Max 

momentum 3-form T a ~ maxenergy3(— a) are defined by 

G a :=±ri a f h AR*' (29) 

and 

Max X 

T a = e Q jL Ma x + (e a \F) AH, with L Max := -- F A H , (30) 
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respectively. The excitation H ~ excit2 and the electromagnetic field strength 
F ~ f arad2 are given by 

H=-H Adt + V, F = E A dt + B (31) 



(Tt magnetic and D electric excitation, E electric and B magnetic field). The 
ubiquitous r?-basis is dual to the i?-basis, that is, see [16], 

Va = X , V*P = * (#c A tip) , etc. (32) 



In a Riemannian space, we can decompose the Einstein 3-form into 2 irre- 
ducible piecesQ: 



G a — - r/apj A 



with 



R:=e a \ep\K#. 



(34) 



If we substitute this into the Einstein equation in (28), then it splits into two 
independent pieces, a tracefree one and the trace: 



1 / . a Max 

-Va^A ®Rfr = 2 T a . 



R = 4X. 



(35) 



Max Max 

Note that T a is tracefree: -d a A T a = 0. It is now evident of how we can directly 
use our decomposition program of the last section in order to verify the field 
equations (35). But, of course, we could also compute the Einstein 3-form 
directly according to its definition (29). 



3.1 P&D coframe o(a) and electromagnetic potential potl 



The P&D solution is described in terms of coordinates x M = (r, q, p, a) ~ 
(tau, q, p, sigma). We use Eq.(3.30) of Plebanski and Demianski [27], see also 
[9,4,5,10]. Then the orthonormal coframe reads, 

1 In a metrioaffine space, the definition (29) kills the strain curvature and also 3 
of the rotational curvature pieces. We find: 

1 6 

G a = -7 ]af3l Aj2 (I) W^. (33) 

1 1=4 

Consequently, in a metric-affine space, besides ricsymf and scalar, only ricanti 
surfaces as an additional quantity in the Einstein 3-form. 
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if 



1 9l 

H\ A 



(dr — p 2 daj , 



1 A 



dp, 



^ = ^^(dT + q 2 da) , 
with the metric 



(36) 



(37) 



or 

9 = JP {"I {dT ~ p2rf °" )2 + I ^ + I ^ + I ( *" + g2 ^ )2 } • (38) 

The functions V, Q, A, and H are polynomial in p and q. With the parameters 
m, n, e Q , g , 6, e, /x, they read: 



7? := (b - g 2 ) + 2np - ep 2 + 2m/ip 3 - /i 2 (ft + e 2 ) + ^ 



Q := (6 + e 2 ) - 2mg + eq 2 - 2n/ig 3 - ^ 2 (& - <?o) + 



P 



A:=p 2 + g 2 , 
H :—l — npq . 



(39) 



We recover the original version of P&D [27], if we put b := 7 — A/6 and /i = 1. 
The parameter //, see [10], generalizes the P&D solution slightly. 



The electromagnetic potential A appropriate for this solution can be expressed 
as follows, see [27] (e Q = electric and g Q = magnetic charge): 



1 



A = X K e ° 9 + 9o p) dr + (g D q-e Q p)pq da] 



A 
H 



e a q ^ + 9oP^ 



A \VQ 



(40) 



In the following we will show of how one can put this metric into Excalc. 
In an interactive Reduce 3.6 session, the Excalc package can be loaded by 
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the command loacLpackage excalc$ . In the next step, Excalc has to be 
acquainted with the rank of the differential forms to be used in the code. This 
is achieved by the pf orm declaration. Subsequently it has to be specified by 
means of f domain on which variables these forms do depend. Afterwards one 
defines the coframe: 

% Specifying the coframe o(a) * 

% file coframe. exi, 1998-03-15 
% no prior input 

loacLpackage excalc$ 

pf orm {hh , sqrtqq , sqrtpp .delta , polynomq , polynomp}=0$ 

f domain hh=hh(p,q) , sqrt qq= sqrtqq (q) , sqrtpp=sqrtpp(p) ,delta= 

delta(p,q) ,polynomq=polynomq(q) ,polynomp=polynomp(p)$ 

coframe o(0) = ( (1/hh) *sqrtqq/sqrt (delta) )*(d tau-p**2*d sigma) , 
o(l) = (( 1/hh) *sqrt (delta) /sqrtqq)* d q, 
o(2) = (( 1/hh) *sqrt (delta) /sqrtpp)* d p, 
o(3) = (( 1/hh) *sqrtpp/sqrt (delta) )*(d tau+q**2*d sigma) 

with 

metric g = -o(0)*o(0) + o(l)*o(l) + o(2)*o(2) + o(3)*o(3)$ 
frame e$ 
$end$ 

Here the frame e a ~ e(— a) is dual to the coframe. 
In accordance with (39), we introduced the functions: 

hh ~ 1 — fxpq, 
delta ~ p 2 + q 2 , 
sqrtqq ~ v^Q, 

sqrtpp ~ VV . (41) 

Here Q and V, see (39), are polynomials of quartic order in q and p, respec- 
tively. For the time being, we will only read in A and H, since they are not 
very complicated: 

delta := p**2 + q**2; 
hh := 1 - mu*p*q; 

However, the explicit form of and \fV will be left open. 
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It is an old trick to treat square roots of functions in the way we did it. Bringing 
in square roots explicitly, will usually slow down a program appreciably or 
even may blow it up beyond the capacity of the computer. We can express the 
derivatives of \fQ and \J~V by means of the chain rule, also their squares get 
a new name: 

@(sqrtpp,p) := @(polynomp,p)/(2*sqrtpp) ; 
@(sqrtqq,q) := @ (polynomq, q) / (2*sqrtqq) ; 

sqrtpp**2 := polynomp; 
sqrtqq**2 := polynomq; 



3.2 The curvature riem2(a,b) of the P&D solution 



Let us first bring in the metric volume element and the 77-basis: 

% The metric volume element eta4 and the eta basis * 
°/„ file eta.exi, 1998-03-15 

°/„ prior input: coframe o(a), frame e(a), metric g(a,b) 

pform eta0(a,b,c ) d)=0,etal(a,b,c)=l,eta2(a,b)=2,eta3(a)=3, eta4=4$ 



= # 1$ 

= e(a) _| eta4$ 

= e(b) _| eta3(a)$ 

= e(c) _| eta2(a,b)$ 

= e(d) _| etal(a,b,c)$ 



eta4 
eta3(a) 
eta2(a,b) 
etal(a,b,c) 
eta0(a,b,c,d) 
$end$ 



This is simple geometry. Now we compute connection and curvature: 

°/ Riemannian connection chrisl(a,b) and curvature riem2(a,b) * 

% file riemann.exi, 1998-03-15 
°/„ prior input; o(a), e(a), g(a,b) 

factor ~,o(0) ,o(l) ,o(2) ,o(3)$ 

Riemannconx chrisl$ 

°/„ because of our conventions, we have to take the transp. of it 
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chrisl(-a,b) :=chrisl(b,-a)$ 

pform -f_riem2(a,b) , curv2(a,b)}=2, -[einstein3(a) , cosmolog3(a)}=3, 
grscalar=0$ 

% antisymmetric riem2; since we want to turn non-Riem. later 

riem2(-a,b) := d chrisl(-a,b) - chrisl(-a,c) " chrisl(-c,b)$ 
grscalar := e(-a) _| (e(-b) _| riem2(a,b))$ 

einstein3(a) := (1/2) * etal(a,b,c) ~ riem2(-b,-c) ; 
cosmolog3(a) := lam * eta3(a)$ 

curv2(a,b) := riem2(a,b)$ 
$end$ 



3.3 Starting an interactive Reduce session 

Let us start such a session by calling Reduce and reading in the co frame and 
the 77 files. Then we specify A, H, etc. and call the Riemann file: 

reduce 

in "cof rame . exi"$ 

in "eta.exi"$ 

delta := p**2 + q**2; 

hh := 1 - mu*p*q; 

@(sqrtpp,p) := @(polynomp,p)/(2*sqrtpp) ; 
@(sqrtqq,q) := @ (polynomq, q) / (2*sqrtqq) ; 

sqrtpp**2 := polynomp; 
sqrtqq**2 := polynomq; 

in "riemann. exi"$ 

We hasten to point out that we copied the commands listed above from the 
original tex-file of this paper directly into a separate Reduce window. Thereby 
our Reduce-Excalc code is on line computer checked and, thus, absolutely 
trustworthy. 

After a few seconds (with a Sun sparcstation) the curvature will be in the 
output in terms of the up to now unspecified functions sqrtqq, sqrtpp, 
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polynomq, polynomp. Nevertheless, if we read in the rotational and the strain 
curvature decomposition programs of above, 

in "rotcurvdecomp . exi"$ 
in "straincurvdecomp . exi"$ 

we immediately find weyl2, ricsymf2, scalar 2 as non- vanishing, whereas 
the whole rest vanishes, as it must be in a Riemannian spacetime. Before 
we look into more details of the surviving curvature pieces, we discuss the 
coupling to the electromagnetic field. 



3.4 The electromagnetic field as source of gravity: maxenergy3(a) 

Our Maxwell program has been published earlier, see [30]. It has to be pre- 
ceded by the coframe statement and the specification of the frame e(a). This 
has been already executed in the last subsection. Therefore we can take the 
electromagnetic potential from (40) and type it in directly into the Maxwell 
program: 

°/„ Maxwell equations metric free * 

'/, file maxwell . exi , 1998-03-15 

°/„ prior input: coframe o(a), frame e(a), metric g(a,b) 
% potential is prescribed 

% field strength farad2, excitation excit2, and the left hand 
% sides of the Maxwell equations are defined 

pform potl=l, {f arad2,excit2}=2, {maxhom3,maxinh3}=3$ 

potl := (l/delta)*( (ee*q+gg*p) * d tau 

+ (gg*q-ee*p) *p*q * d sigma) ; 

farad2 := d potl; 
maxhom3 := d farad2; 
excit2 := # farad2; 
maxinh3 := d excit2; 

°/„ Maxwell Lagrangian and energy-momentum current are assigned 

pform lmax4=4, maxenergy3(a)=3$ 

lmax4 := -(1/2) * farad2 " excit2; 
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maxenergy3(-a) := e(-a) _| lmax4 + (e(-a) _| farad2) " excit2; 
$end$ 

In continuing our interactive Reduce session, we read in this file: 
in "maxwell . exi"$ 

Both Maxwell equations are fulfilled. The electromagnetic energy-momentum 
is determined. If we type in 

potl ; 

the frame version of A will be displayed, cf. (40). The Einstein equation with 
cosmological constant and electromagenetic field source can be written down: 

pform einsteinmax3(a)=3$ 

einsteinmax3(a) :=einstein3(a) + cosmolog3(a) - 2*maxenergy3(a) ; 

The left hand side, namely einsteinmax3(a) if evaluated, has already a very 
compact form. Then, substituting the polynomials V and Q of (39), 

7„ substitution 

polynomp := bb - gg**2 + 2*n*p - epsi*p**2 + 2*mu*m*p**3 
-(lam/3 + mu**2*(bb + ee**2) )*p**4; 

polynomq := bb + ee**2 - 2*m*q + epsi*q**2 - 2*n*mu*q**3 
-(lam/3 + mu**2*(bb - gg**2) )*q**4; 

einsteinmax3(a) := einsteinmax3(a) ; 
showtime ; 

the Einstein equation turns out to be fulfilled, too. 

3.5 Conformal curvature two-form weyl2(a,b) of the P&D solution 

Consequently the P&D metric is a solution of the Einstein-Maxwell equa- 
tions, indeed. The P&D gravitational field is characterized by the conformal 
curvature 2-form 

C al3 := (1) R a(3 = l - C^p fAf. (42) 

Therefore it is desirable to put this in a compact form. Our decomposition 
program ran already. Thus we only have to collect the result by turning the 
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expansion off and calling for the Weyl 2-form: 
off exp$ 

weyl2(a,b) :=weyl2(a,b) ; 
We find thereby (Heinicke [18]): 

C aP = f(a,/3) {-A$ a A^ + B*(f A^)} , (43) 



with 



A . = (wq - i 



B:- 



(p2 + q 2f 

(fipq - l) 3 



(p 2 + Q 2 Y 
and 



m (3p 2 — q 2 ^j q + n (} 


> 2 - 3g 2 ) 


m (3q 2 — p 2 ^j p — n (i 


I 2 - 3P 2 ) 



f ^ r) 2 for (a, 0)e {(0,1), (1,0), (2, 3), (3, 2)} 

-1 else 



Of course, the Weyl 2-form is a real quantity. Accordingly, we can also repre- 
sent it by its self-dual part, see [20] (i— imaginary unit), 

+ C a(3 := l - (C a/3 + i *C a p) = l - + C, vaP fAf, (45) 



or, in computer code, by 
on complex$ 

pform self dualweyl2(a,b)=2$ 

selfdualweyl2(a,b) := (1/2) * (weyl2(a,b) + i * #weyl2(a,b)) ; 

It can be very compactly displayed [18], 

+C a(B = f(a, (3) C {#* A ^ + i* (iT A/)} , (46) 



where 



C := --{A + iB) = - — — T \ n - im - e + g ) — — . (47) 



2 V ' 2{p-iqY \ V u ™> p + iq 
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Often it is thought that such a compact representation is only possible in a 
null frame (i.e., by means of a Newman- Penrose tetrad). Obviously, however, 
this is not true. It is an effect of the power of the moving frame aspect of 
exterior calculus, not of that of a null coframe. 

By the same token, we computed also the two quadratic curvature invariants 
*(W a p A* W af3 ) and *(W a p A W af3 ), but their discussion is beyond the scope 
of this paper. 



4 Torsion and nonmetricity and their irreducible decompositions 



In a metric-affine spacetime, the field strengths, besides the curvature two- 
form curv2(a,b), are given by the components of the torsion two-form 

torsion2(a) ~ T a = D$ a = d$ a + V p a A ^ = ]-T^ a 0" A , (48) 



and the nonmetricity one-form 

nonmetl(-a, -b) ~ Q afS = -Dg afS = -dg afS + T^g^ + 1^ 1 g ai 

= Q ia( sdx l . (49) 



Although torsion and nonmetricity are genuine field strengths, they can be 
reinterpreted as parts of the connection. The linear connection can be ex- 
pressed in term of metric, coframe, torsion, and nonmetricity. 



4-1 Torsion torsion2(a) and its three irreducible pieces 

For the torsion, we have the following irreducible decomposition 

T a = (i) T a + (2) T « + ( 3 ) T « 

~ tentor + trator + axitor , (50) 



where 



trator2(a) ~ —i} a A (ep\T^j , (51) 
axitor2(a) ~ - e a \ (dp A T' 3 ) , (52) 
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tentor2(a) ~ T a - (2) T° - (3) T a , 



(53) 



or, in other 'words', 

% Irreducible decomposition of the torsion torsion2(a) * 

°/„ file torsiondecomp . exi , 1998-03-15 

% input: o(a), e(a), g(a,b) , torsion2(a) 

pf orm {tentor2 (a) , trator2 (a) , axitor2 (a) }=2$ 

trator2(a):= (l/3)*o(a) ~ (e(-b) _| torsion2(b) ) ; 

axitor2(a):= (l/3)*e(a) _| (o(-b) ~ torsion2(b) ) ; 
tentor2 (a) : = torsion2 (a) -trator2 (a) -axitor2 (a) ; 
$end$ 



4-2 Nonmetricity nonmetl (a, b) and its four irreducible pieces 

The Weyl one-form Q := Q a a /4: = —g al3 Dg a 0/4: is one irreducible piece of the 
nonmetricity. Thus we find the traceless part of the nonmetricity as 



% a p '■= Qap - Qg a p- (54) 

The traceless nonmetricity (its 'deviator', in the language of continuum me- 
chanics) can be further 'contracted': 

A Q :=e^ Q/3 , A := A a $ a , (55) 

e Q :=*(# Q/3 A^), 8:=fAe a , Q a := Q a - ^e a \Q . (56) 

From these quantities we can build up the four irreducible pieces: 

binoml(-a, -b) ~ - * ($ {a A , (57) 

4 / 1 \ 

vecnoml(-a, -b) ~ - h3 w A a) - -g af3 A) , (58) 

conoml(-a, -b) ~ g a pQ, (59) 

trinoml(-a, -b) ~ Q aP - ^Q ap - ^Q^ - ^Q ap . (60) 



19 



trinom represents a tensor of 3rd rank and binom one of 2nd rank, with vecnom 
and conom we denote the covector and the vector pieces, respectively. Clearly 
conom is equivalent to the Weyl covector. They add up according to 



Qa/3 — (l> Qal3 + + + W Qa/3 

~ trinom + binom + vecnom + conom . 



(61) 



% Irreducible decomposition of the nonmetricity nonmetl(a,b) * 

% file nonmetdecomp . exi , 1998-03-15 

% input: o(a), e(a), g(a,b) , nonmetl(a,b) 

pform lamzero(a)=0, {weylcovectorl , nomtracef reel (a, b) , lamonel , 
binoml(a,b) ,vecnoml(a,b) ,trinoml (a,b) , conoml (a,b) }=1 , 
{thetatwo2 (a) , omega2 (a) }=2 , thetathree3=3$ 



weylcovectorl 

nomtracef reel (-a, -b) 

lamzero(-a) 

lamonel 

thetatwo2(-a) 

thetathree3 

omega2(-a) 



= nonmetl(-c,c)/4$ 

= nonmetl (-a, -b) - g(-a, -b) *weylcovectorl$ 

= e(b)_ | nomtracef reel (-a, -b)$ 

= lamzero(-a)*o(a)$ 

= #(nomtracef reel (-a, -b) ~o (b) ) $ 

= o(a)~thetatwo2(-a)$ 

= thetatwo2 (-a) - (1/3) *e (-a) _ I thetathree3$ 



binoml(-a,-b) := (1/3) *#(o (-a) ~omega2(-b)+o (-b) "omega2 (-a) ) ; 
vecnoml(-a,-b) := (4/9)*((o(-b)*lamzero(-a)+o(-a)*lamzero(-b))/2 

-g(-a,-b)*lamonel/4) ; 
conoml (-a, -b) := g(-a, -b) *weylcovectorl ; 
trinoml(-a,-b) := nonmetl (-a, -b) - binoml (-a, -b) 

-vecnoml (-a, -b) - conoml (-a, -b) ; 

clear lamzero(a) ,nomtracef reel (a,b) , lamonel ,thetatwo2 (a) , 
omega2(a), thetathree3=3$ 

$end$ 



4-3 The post-Riemannian covector triplet 



With propagating nonmetricity Q a p two types of charge are expected to arise: 
One dilation charge related via the Noether procedure to the trace of the non- 
metricity, the Weyl covector Q = Qidx\ It represents the connection associ- 
ated with gauging the scale transformations (instead of the £/(l)-connection 
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in the case of Maxwell's field). Furthermore, nine shear charges are expected 
that are related to the remaining traceless piece ($ a p := Q a p — Q g a p of the 
nonmetricity with 4 x (4 + 4+1) components. 

For the torsion and nonmetricity field configurations, we concentrate on the 
simplest non-trivial case with shear. According to its irreducible decomposi- 
tion (61), the nonmetricity contains two covector type pieces, namely ^Q a p = 
Q g a/ 3, the dilation piece, and 

{3) Q a p = \ ($ {a e 0) \ A - , with A = 0Vj# a/9 , (62) 

a proper shear piece. Accordingly, our ansatz for the nonmetricity reads 

Q a( s = {3) Q a( s + {4) Q a( s ■ (63) 

The torsion, in addition to its tensor piece, encompasses a covector and an 
axial covector piece. Let us choose only the covector piece as non-vanishing: 

r = (2) T a = -rAT, with T:=e a \T a . (64) 

o 

Thus we are left with the three non-trivial one-forms Q, A, and T. We shall 
assume that this triplet of one-forms shares the spacetime symmetries, that 
is, its members are proportional to each other. 



5 Applying the decomposition programs in metric-afRne gravity: 
The P&D metric together with a non-Riemannian covector triplet 
as exact solution of a dilation-shear Lagrangian 

In order to explore the potentialities of metric-afline gravity, see [11], and 
of the corresponding computer algebra programs, we will choose the simple 
non-trivial dilation-shear gravitational Lagrangian, 

Km-sh = ~ (R a(B A Vaf) + PQ A *Q + 7 T A *T) - R a a A * V > ( 65 ) 

with the dimensionless coupling constants a, (3, and 7. We will choose units 
such that l 2 — 1. We add the Maxwell Lagrangian, see (30): 

L = Vdil-sh + -^Max • (66) 



21 



Explicitly we have: 



L = ~ w A Vaf} + p Q A * g + 7 T A * T ) " I a Rc * a A * Rf)fi 

— i F A*F . (67) 

Note the formal similarities of the curvature square Lagrangian with the seg- 
mental curvature \R a a = dQ and the Maxwell Lagrangian with the electro- 
magnetic field strength F = d A. 

We will search for exact solutions of the field equations belonging to this La- 
grangian. Since the matter part Lmsx does not depend on the connection IV 3 , 
the hypermomentum A"^ := 5LM ax /^r Q , /3 vanishes, A a p = 0, and the only 

Max 

external current is the electromagnetic energy-momentum current S a —T a . 
We will only be able to find non-trivial solutions, if the coupling constants 
fulfill the constraint 

7 = -?-?_. (68) 



This is the best we can do so far. 



It has been shown by Obukhov et al. [24] - - for the related work of Tucker 
et al., see [40] - how the general solution belonging to Lagrangians of the 
type (67) can be constructed: One starts with an electrovac metric of general 
relativity and the corresponding electromagnetic potential. We will take the 
P&D metric (38) and its potential (40). Then the post-Riemannian triplet of 
(63,64) is constructed which is patterned after the electromagnetic potential, 

Q = k u , A — kiu, T = k 2 u , (69) 



with the 1-form 



u=^[(N c q + N g p) dr+ (N g q- N c p)pqd(T} 

H f^ q ^ + MA. (70 ) 



In Excalc this input reads: 

% Triplet ansatz for nonmetricity and torsion * 
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% file triplet. exi, 1998-03-15 

°/„ prior input coframe o(a), frame e(a), metric g(a,b) 



pform {omegal, qql, llaml, ttl}=l$ 

omegal := (l/delta)*( (Ne*q+Ng*p) * d tau 

+ (Ng*q-Ne*p) *p*q * d sigma) ; 



qql 

llaml 

ttl 



= kO * omegal$ 
= kl * omegal$ 
= k2 * omegal$ 



clear omegal$ 
$end$ 

These are the three one-forms which are fundamental for the exact solution. 
Now these one-forms have to be assigned to the corresponding pieces of the 
torsion T a ~ torsion2(a), the contortion K a p ~ contorl(— a, — b), and the 
nonmetricity Q alS ~ nonmetl(a, b): 

°/„ Post-Riemannian assignments * 

% file postriem.exi, 1998-03-15 

°/„ prior input o(a), e(a), g(a,b,), qql, llaml, ttl 

pform torsion2(a)=2,{contorl(a,b) ,nonmetl(a,b)}=l$ 

torsion2(a) := (1/3) * o(a) ~ ttl$ 

nonmetl(a,b) := (2/9)*( o(b) * (e(a) _| llaml) 

+o(a) * (e(b) _| llaml) 
-(l/2)*g(a,b)*llaml ) + qql*g(a,b)$ 

contorl(-a,-b) := (l/2)*(e(-a) _| torsion2(-b) 

e(-b) _| torsion2(-a) ) 
-(l/2)*(e(-a) _| (e(-b)_|torsion2(-c)))"o(c)$ 

$end$ 

The constants k , k±, k 2 in front of the triplet can be expressed in terms of the 
coupling constants of the dilation-shear Lagrangian: 

kn = —— , ki = — — — — , ko = 6 . (71) 
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N c and N s are the quasi-electric and the quasi-magnetic dilation-shear-spin 
charges, respectively. Now, in metric-affine gravity, the polynomials depend 
also on these quasi-charges in a fairly trivial, i.e. in an additive way, 



Q 

A 
H 



= {b- gl - Gl) + 2np - ep 2 + 2m/ip 3 - [/i 2 (b + e 2 + £ 2 ) 
= (b + e 2 + Si) - 2mq + eq 2 - 2n M 3 - [/i 2 (b - g 2 Q - Gl 



2 , 2 

=p +q , 

= 1-Hpq. 



(72) 



The post-Riemannian charges ^ c and N g are related to the post-Riemannian 
pieces S Q and Go, entering the polynomials, according to 



/oi fa 
-N e , Go = k ^j-N g , with a>0. 



(73) 



The field equations of metric-affine gravity can be expressed in term of the 
irreducible components of curvature, torsion, and nonmetricity - and these 
quantities we have already computed. With a corresponding computer pro- 
gram we verified the correctness of the exact solution discussed in this sec- 
tion. The presentation of these programs, however, we will leave to another 
occasion. 



6 Discussion 



Effective work in gravity can nowadays only be done if one has the standard 
computer algebra methods and packages available. We hope to have shown 
that the Reduce-Excalc tool is very efficient in the context of the Einstein- 
Maxwell equations as well as in gauge theories of gravity, like metric-affine 
gravity, where the spacetime is a non-Riemannian manifold. 
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